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Abstract. We consider how measurements of the 21cm radiation from the epoch of reion- 
ization (z = 8 — 12) can constrain the amplitudes of various 'shapes' of primordial non- 
Gaussianity. The limits on these shapes, each parametrized by the non-linear parameter /nLj 
can reveal whether the physics of inflation is more complex than the standard single-field, 
slow-roll scenario. In this work, we quantify the effects of the well-known local, equilateral, 
orthogonal and folded types of non-Gaussianities on the 21cm power spectrum, which is ex- 
pected to be measured by upcoming radio arrays such as the Square-Kilometre Array (SKA). 
We also assess the prospects of the SKA in constraining these non-Gaussianities, and found 
constraints that are comparable with those from cosmic-microwave-background experiments 
such as Planck. We show that the limits on various /nl can be tightened to 0(1) using a 
radio array with a futuristic but realistic set of specifications. 
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1 Introduction 

The statistical nature of primordial density perturbations can reveal much about the physics 
of the early Universe. If the simplest theory of single-field slow-roll inflation is correct, the 
perturbations should obey an essentially Gaussian statistics [1]. Any detection of primordial 
non-Gaussianity would hint at novel physics during inflation, or perhaps even an alternative 
theory of the early Universe (see [2, 3] for reviews). 

Much effort in the search for primordial non-Gaussianity has been directed towards 
constraining the 'local' non-Gaussianity parameter, /nl, which occurs in the expansion of 
the nonlinear Newtonian potential, as a series of Gaussian random field, (j), 

$ = </' + /nl('/'' -(</>')) + ... (1.1) 

The most stringent constraint on /nl comes from the cosmic microwave background (CMB) 
anisotropics which give /nl = 32 it 42 {2a) [4]. 

In addition, there may be other kinds of non-Gaussianity, including the so-called equi- 
lateral, orthogonal and folded non-Gaussianities which, as we explain in more detail in the 
next section, are signatures of specific physics in the inflationary era. We refer to these va- 
rieties as shapes of non-Gaussianity, as they are associated with different 3-point correlation 
functions, i?(ki, k2, ks), which peak when the Fourier wavevector, kj, of the perturbations 
form different triangular configurations [5-7]. Current constraints on these shapes are much 
weaker than that of the local model. At 2(7, the constraints from WMAP are /^x^'' = 26ib240, 
/nl*^ = —202 lb 208 [4] (there is no published constraint on /^l^'^'^ present). These con- 
straints are expected to be improved by at least an order of magnitude with the highly 
anticipated results from the Planck satellite^, as well as constraints from various upcoming 
large-scale-structure (LSS) surveys (see e.g. [8] for a review). 

More recently, it has been shown that non-Gaussianity could leave a detectable signature 
in the 21cm radiation due to spin- flip transitions of neutral hydrogen during the epoch of 
reionization or earlier during the cosmic Dark Ages [9-14]. In particular, Joudaki et a/. [10] 

^www. rssd. esa. int/planck 
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showed that there are good prospects of constraining local /nl with the Square-Kilometer 
Array (SKA)^ to within |/nl| ^ 10 using measurements of the 21cm power spectrum from 
redshift z ~ 7. In this work, we extend this pursuit to other shapes of non-Gaussianity. 
In particular, we investigate the prospects of the SKA in constraining /nl associated with 
different non-Gaussianity shapes, in comparison with the CMB or LSS constraints. 

Throughout this work we assume in our fiducial model the following cosmological param- 
eters [15]: present density parameter for matter(r2m = 0.276), baryon = 0.046), radiation 
(Qr = 5x 10-^), dark energy (Ode = 1 — — f^r); scalar spectral index, iig = 0.96, and 
normalization of matter power spectrum, erg = 0.81. 

2 Shapes of non-Gaussianity 

The bispectrum, B^, is defined as the three-point correlation function of the Newtonian 
potential, <I>. In this work, we shall consider four common forms of the bispectrum described 
below. 

1. The local model is the dominant type of non-Gaussianity in single field inflation and 
curvaton models [16-18]. The bispectrum is given by 

i?'-(ki,k2,k3) = 2/;°£^'-, (2.1) 
jc-ioc ^ p^p^ _^ 2 perm. (2.2) 

where Pi = P^{ki) oc k^""^ is the primordial power spectrum of In wavevector space, the 
bispectrum peaks when the vectors ki,k2,k3 form a 'squeezed' triangle in which /ci ~ A;2 ^» 
k3. 

2. The equilateral model of non-Gaussianity is dominant in inflation models with higher- 
derivative terms in the Lagrangian. The bispectrum peaks when ki, k2, ka form an equilateral 
triangle [5]. The precise expression for the bispectrum is model-dependent and may take a 
complicated form, but for computational purposes they can be approximated by a 'template' 



which comprises factorisable pieces of the power spectrum: 

B'^^'ikiM, ka) = 6f^t {-F'°' - + T"") , (2.3) 

= {PiP2P3f\ (2.4) 

-pB ^ p^/^p^/^p^ + 5 perm. (2.5) 



(see [19] for numerical comparison between the templates and physical shapes.) 

3. The orthogonal model was constructed to be orthogonal to the equilateral shape with 
respect to an inner product of bispectra as deflned in [7]. It also occurs in models with 
higher derivatives in the Lagrangian. The template for the bispectrum is 

B^'-^^ki, k2, kg) = 6/^2^ (-3-F'°^ - 8J-^ + 3-F^) , (2.6) 

which peaks on equilateral and folded conflgurations (see below). 

"www . skatelescope . org 
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4. The folded model of non-Gaussianity is prevalent in inflation models in which the initial 
Bunch-Davies vaccum state is modified [6]. The template for the bispectrum is 

i?f°l(ki, k2, kg) = ^ [b''^' - , (2.7) 

peaking on a flattened triangle configuration with ki = 2k2 = 2/C3. 
3 A model of reionization 

We are looking for signatures of the various non-Gaussianity shapes in the 21cm signal from 
the epoch of reionization (z ~ 8 — 12). The intensity of the 21cm radiation depends on the 
local ionized fraction, Xj. We define the ionized- fraction contrast as 

4 = (3.1) 

Xi 

where Xi is the mean ionized fraction. Both Xi and Xi depend on the model of reionization. 

In this work, we appeal to the reionization model of Alvarez et al. [20], in which the 
ionized fraction Xi is a function of the local overdensity, 5, and redshift, and satisfies 

ln[l - x,{6, z)] = -C(5, ^)/coii(5, z, M), (3.2) 

where C is the reionization efficiency and /coU is the collapse fraction in a region associ- 
ated with a top-hat smoothing function with smoothing scale R (or smoothing mass M = 
AttR^Po/S). Alvarez et al. noted that this simple parametrization reproduces the behaviour 
expected when the number of ionizing photons is large (C/coii ^ 1) and small compared with 
atoms (see their Appendix for detail). A particularly simple form of C can be obtained in 
the "Stromgren" limit in which a photon is emitted for every HI atom formed upon recom- 
bination. In this case, it can be shown that 

m = -^^^Co{l-6), (3.3) 

with Co ^ constant and \5\ <^ 1, which holds in our application. 

The collapse fraction depends on the halo formation model. For example, in the ex- 
tended Press-Schechter theory [21, 22], /coU depends on the difference between the local 
overdensity 6{z) and the critical collapse overdensity 6c{z)^. 

where Mmin is the minimum mass for collapse, which we take to correspond to the virial 
temperature of 10^ K (at which point atomic hydrogen line cooling is efficient) 

M^i„ = 2.77xlO^/."iMe(^) (^) . (3.5) 



^We shall make a strict distinction between 5c and 5c(z): Sc{z) = 5c/D{z) and 5{z) = 5/D(z), where 5c 
and 5 are overdensities extrapolated to z = and D{z) is the growth function. 
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Following Tashiro & Ho, we approximate the relation between the matter density con- 
trast, 5, and the ionized- fraction contrast 6x by 



6x = bxS, (3.6) 

where bx is referred to as the bias parameter for the ionized fraction. We can obtain a relation 
between bx and the usual halo bias, b, by expanding (3.4) about small \6\ and equating the 
result with /coU = /coii(0)(l + {b - 1)6), yielding 

bx = ^^Hl-Xi){b-2), (3.7) 
where b is the average bias of dark matter halos with mass above m^in 

6=1 + W — — , (3.8 

with z/ = 6c{z)/a and z^min = 6c{z) / amin- The expression for the mean ionized fraction can 
be obtained by setting (5 = in (3.2), 



ln(l - Xi) = -Co erfc (^t'min/v2j . (3.9) 

In summary, the ionized- fraction bias, b^, can be calculated once a history of reionization is 
known. 

As pointed out in [20], the average halo bias b can also be calculated by averaging the 
mass-dependent Eulerian bias over the multiplicity function f{u) (see [23] for details). 

fps{i^)bpsii^)dlniy 
/-min^Ps(^)rfln^ 

^ps = l + ^— , (3.11) 

6c 

/ps(z^) = i^e-""/'. (3.12) 

This gives us a way to extend the bias calculation to other collapse formalisms. We 
outline the calculation using the Sheth-Tormen formalism [24] in the Appendix, and comment 
on the results in section 5 

In the next section we shall see how the various non-Gaussianity shapes manifest in the 
ionized-fraction bias. 

4 The 21cm power spectrum in the presence of non-Gaussianities 

The intensity of the observed 21cm signal can be measured by the differential brightness 
temperature, 6Th (the difference between the 21cm and CMB temperatures). In the optically 
thin limit, it can be shown that the mean brightness temperature is given by (see e.g. [25]) 

AT^r^-o-^ss T^- ^ Tcmb\ n^h^ f 0.15 1 + zV/^ 

5r,(z)« 23.88 mKxHi(^l-^J^(^^^J , (4.1) 

where xhi is the average neutral fraction of hydrogen (rrni + Xi = 1). 
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Assuming the fluctuations are linear, one can show that the 21cm power spectrum can 
be expanded as [26, 27] 



(4.2) 



where /i is the cosine of the angle between wavevector k and the line-of-sight vector n. Here 



we define 



6TjjX^iPL, where -Pl(A;) is the linear power spectrum oc M.'^{k)kP"' ^, and 



M{k) 



2k'^T{k)D{z 



(4.3) 



where T{k) is the transfer function and D{z) the growth factor. The remaining components 
of (4.2) are defined by 



^ x5 — ^x ^ i5(5 ) 
XX — ^ (SiS 1 



(4.4) 
(4.5) 



which are the density-ionization cross-correlation and the ionization autocorrelation respec- 
tively. These expressions reduce the 21cm power spectrum to 



i^5T(k) 



2\2 



[1 - fex + [l') 



(4.6) 



Non-Gaussianity affects the 21cm power spectrum by introducing a characteristic scale de- 
pendence in the bias hx- This is analogous to the non-Gaussian effects on the halo bias, 
which has already allowed a number of cluster surveys to place limits on non-Gaussianity, 
particularly the local type (see e.g. [8, 28] and references therein). Calculations of the halo 
bias in the presence of other shapes of non-Gaussianity were explored by Matsubara [29] 
(amongst others [8, 19, 30]). The calculations in this work are based on the analytic results 
of [29]. 

Writing 6^^ = ^Gaussian _^ Matsubara showed using his theory of "Integrated Per- 
turbation Theory" that in the large-scale limit, the bias shift due to non-Gaussianity in the 
PS formalism is given by 



A6ps(fc,M) = i<5,(6ps - 1)X(A:,M) + \-^^. 

I Z am aM 



(4.7) 



The function X(/c, M) is derived from the bispectrum shapes as follows 

1 



I{k, M) 



4/nl 
M{k) 



X < 



3A;" 
3/^71 /2 
-3A:7i 



(1 + 05)^1 



(local) 
(equilateral) 
(folded) 
(orthogonal) 



(4.8) 



where 



as 
lt{M) 



3 
1 



(4.9) 



dk_ 
1 



k^-W\kR)PL{k). 
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Here W{kR) is the Fourier transform of the top-hat window function associated with mass 
M = 4TrpoR'^/3. 

The main idea is to replace the bias b in the Alvarez bias (3.7) by the non-Gaussian 
version, which we calculate as follows. First, we substitute the Gaussian bias by the non- 
Gaussian expression 

bps = bps + A6ps(A:,m), (4.10) 

where A6ps is given by the Matsubara expression (4.7). This allows us to calculate the 
mass-averaged bias by an averaging analogous to (3.10): 



TOO dn jj^ 



where we have transformed the dv integral to an equivalent one in dM. The mass function, 
dHlU' satisfies 

dn _ pm ^ dlng-^ _ f{v) dv 
dlnM ~ M^^""^ dlnM " ^" u dM' ^ ^ 

However, the mass function is also affected by the presence of non-Gaussianity. Here, we 
modify the mass function using the prescription of [31], in which the Gaussian mass function 
is multiplied by a correction factor T^-ng^ 



dn dn 
dlnM " ^^dlnM 



7^NG = exp 



Gaussian 

+ A 



6l_ dSs 
6 A din (7 A/ 



(4.13) 
(4.14) 



A^yi-^, (4.15) 
where the third cumulant, 53(M), satisfies the relation (see e.g. [8, 32]) 

a^Ss{M) = C kldki [\jdk2 C diiB^{kiMM)T\A{h,R), (4.16) 
JO Jo J~i ~~[ 

ks = {kl + ki + 2fikik2)^^\ A{ki,R)=M{ki)W{kiR). (4.17) 

Putting all these ingredients together, we then obtain the (scale-dependent) non-Gaussian 
bias in the ionized fraction 

6NG(^) = ^i^ln(i_Xi)(6(A;)-2). (4.18) 

Xi 

Figure 1 shows the ionized-fraction bias shift as fraction of the Gaussian bias, 

A !,NG ^Gaussian 

r = 12 f4 19) 

^ ^Gaussian ' \ ■ / 



with Ta; = when /nl = 0. The figure shows the results for various shapes of non-Gaussianity 
with /nl = 50, Co = 50, evaluated at 2; = 8. The curve for the orthogonal type is always 
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10"^ — 
0.01 



0.1 



k/h"''Mpc 



Figure 1. The fractional cliange in the ionized- fraction bias = A6i:/6j. Gaussian due to non- 
Gaussianity of various shapes, calculated at z = 8. The parameters in panel (A) are /nl SOi 
ionization efficiency Co — 50. 

negative and its absolute value is shown here. These curves are reminiscent of those for 
the scale-dependent halo bias seen by previous authors [19, 29, 30]. The strongest scale- 
dependence is exhibited by the local type. The folded type exhibits not only a weak scale- 
dependence of the bias, but also a low amplitude of the bias shift, and thus we can expect 
the poorest constraints for the folded-type non-Gaussianity. 

We note that our approach is different from those in [10, 12], which treated the bias-shift 
in the ionized fraction as an extrapolation of the formula by Dalai et al.[28] 



where ~ 1 is the critical density for collapse in an ionized region. This extension effectively 
assumes that there is a single physical mechanism governing the clustering of linear density 
fluctuations and that of ionized hydrogen clouds, even in the presence of primordial non- 
Gaussianity. This assumption demands further tests from iV— body simulations with non- 
Gaussian initial conditions, which have yet to be performed on such scales, nor has it been 
verified for other types of non-Gaussianity. Our calculations do not rely on this extrapolation, 
and is therefore more self-consistent. 

5 Fisher-matrix analysis 

We now perform a simple Fisher-matrix analysis to assess the potential of upcoming radio 
arrays such as the SKA in constraining non-Gaussianity using the 21cm power spectrum. 

Radio interferometers do not measure the power spectrum directly in k space. Instead, 
they measure the complex visibility function, V, in array-configuration space u = {u,v,u^^), 
where the {u, v) plane is perpendicular to the line of sight and the values of u and v depend 
on the geometry of the array. In the small-angle approximation, we can write (see e.g. [33]) 



Ah 



'x 




(4.20) 
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A^ant 


Lnnn (m) 


Rc (m) 


FoV (deg^)) 


Ae (m^), z = 8/10/12 


SKA 


1400 


10 


330 


232 


50/77/104 


Futuristic 


10^ 


1 


1000 


20626 


1/1/1 



Table 1. Parameters adopted for the Fisher forecasts for the SKA and a futuristic radio array. 

where A is the relative antenna area (equals unity in the line-of-sight direction n). Eq. 5.1 
can be thought of as the Fourier transform of the 21cm differential brightness temperature 
(weighted by A) into the configuration space. 

In the {u,v) plane, we can write u_l = b/A2i(z), where b is a baseline vector between 
two receivers. Hence, u^/27r equals the number of wavelengths that fit between two receivers 
in the u and v directions. The component can be calculated from the shift in frequency 
of the signal from the central redshift of a particular redshift bin. 

The components of u along and perpendicular to the line of sight are related to those 
of the Fourier wavevector k by 

(u±,n||) = {dA{z)\<iuy{z)k\\), (5.2) 

where dA is the angular diameter distance to redshift z and y{z) = 21cm(l + z)'^ /H{z). In 
terms of the baseline L, u± = 2ttL/ X(z), X2i{z) = 21(1 + z) cm. Using these relations, we 
can deduce that the power spectrum in u space is given by 

PstH = :^P5r(k). (5.3) 

The noise power spectrum is given by [34, 35] 

A'^Tsvs \ 1 



where the system temperature Tgys = 280((1 + z)/7.4)^-^ K, to is the observation time, Ag is 
the effective collecting area and n{u±) is the number density of baselines that can observe 
mode u. For surveys with outer core radius Rc and inner nucleus radius Rn, the baseline 
number density can be expressed as 

n{u^) = 27rA^i /"^' rP{r + ^^l^)P{r)dr, (5.5) 
Jo 27r 

'K 0<r < Rn, 

K{RM/rf Rn<r<Rc, (5.6) 



Fir) 



otherwise. 



with normalization constant K = 4e"'^/^A''ant/(5vrii^). 

The error in the measurement in each pixel is given by 



5P(U) = ^^Hu) 



where 

2 



iV=(A'j sineV{z)M/\k, (5.8) 
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Template 



a(/NL) SKA 



ct(/nl) Fut. 



Local 
Equilateral 

Folded 
Orthogonal 



2 

13 
30 
15 



0.1 

1 
3 
1 



Table 2. Forecasts on the /nl constraints of various shapes using the SKA and Futuristic specifica- 
tions in Table 1, and assuming the Press-Schechter collapse formalism. 

is the number of cells in an annulus defined by Ak and spanning half the sphere, and 
V{z) = d\yBx¥6\f , where FoV is the area of the field of view. B (sometimes called the 
bandwidth) is related to the line-of-sight distance, Ar, from the central redshift of a particular 
bin hy B = Ar/y{z). 

Finally, the Fisher matrix is given by 



where the summation extends over the valid pixels in k space. This assumes a Gaussian 
uncertainty in each parameter p (see [36] for an alternative treatment). Assuming further 
that the noise correlation across redshift bins can be neglected, we can add the Fisher matrices 
evaluated at various redshifts. 

In this work, we shall consider two setups, one for the SKA and another for a futuristic 
radio array (such as one considered in [35, 37]). The specifications for these setups are given in 
Table 1. We consider a simple parameter space consisting of the non-Gaussianity amplitude 
and the ionization efficiency, i.e., pa = (/nl^Co)- The background cosmology is fixed using 
the best-fit parameters outlined in the Introduction. Whilst including more parameters in 
the Fisher forecast would most likely yield less stringent constraints on these parameters, 
previous studies taking a similar approach have estimated that when including a full set of 
cosmological parameters, the local /nl constraints would be affected only by ~ 10% [10, 36]. 
Thus we believe that the forecasts we obtain here are indicative of the constraints future 
radio arrays would be able to achieve. 

We proceed to calculate the Fisher matrix as follows. Firstly, a set of experimental 
specifications is chosen (see Table 1) assuming a fixed redshift {z = 8/10/12). We set the 
fiducial values of the parameters to be Co = 50, /nl = and calculate the Fisher matrix over 
the range k G [A^jum) ^max] <md 6 G [^mirn ^max]- The value of femin — 277 /yB is set by the largest 
scale beyond which the signal is dominated by foregrounds (with B = 8.3/7/7.5 MHz), and 
^max = 0.2 Mpc determines the shortest wavelength below which linear calculations break 
down (see [10, 35] for studies of the variation in /cmax)- For a given k, we determine the 
range of 9 such that the mode can be observed by the survey, i.e. those values of 6 such that 
u± = dAksinO. 

The Fisher forecasts on /nl of various shapes are shown in Table 2. Our constraint 
on the local shape is cr(/NL) = [2,0.1] for the [SKA, Futuristic] configurations. The SKA 
constraint in comparable with those expected from Planck. The equilateral and orthogonal 
constraints are similar in both the amplitudes and the scaling and thus distinguishing between 
these shapes using the power spectrum will be challenging. The folded orthognal fared worst 
as expected, with the error roughly twice that of the orthogonal shape. For the latter 3 shapes. 




(5.9) 



-9- 



Template 


a(/NL) SKA 


ct(/nl) Fut. 


Local 


4 


0.3 


Equilateral 


23 


2 


Folded 


56 


5 


Orthogonal 


28 


2 



Table 3. Forecasts on the /nl constraints of various shapes using the Sheth-Tormen formalism. 

the constraints are poorer than the local shape since the region of the strongest signals are 
on large scales which are dominated by foregrounds. In all cases, the SKA constraints will be 
improved by at least an order of magnitude with the Futuristic array, down to (t(/nl) = C'(l) 
for all shapes. 

Our results suggest that there are realistic prospects of contraining the shapes of non- 
Gaussianity with future 21cm power-spectrum measurements. Although the 21cm constraints 
on /nl will not be markedly different from those from upcoming CMB experiments, 21cm 
experiments are complementary to the latter as they constrain /nl on a much smaller scale. 
A multi-scale attack on non-Gaussianity is required if the contributions from the various 
shapes were to be disentangled, and this can be achieved by combining non-Gaussianity 
constraints from the CMB, galaxy clusters and 21cm measurements. Such approach will also 
shed light on the possible scale-dependence of non-Gaussianity [38, 39]. 

We also carried out the Fisher-matrix analysis using the Sheth-Tormen formalism (see 
Appendix A). The values of ^(/nl), shown in Table (3), are roughly twice that of the PS 
formalism. This is in line with expectation from previous studies comparing cluster biases 
in the PS and ST formalisms [14, 29, 40], which found that the amplitude of the bias in the 
ST formalism is significantly less than that in the PS formalism. Alternatively, one could 
also extend our calculations to other functions based on more recent iV-body simulations 
[41, 42], but the bias shift (4.7) will be more difficult to obtain analytically. Nevertheless, we 
expect such calculations to produce /nl constraints that are compatible with the PS and ST 
analyses. 

Our constraints rely on the assumption that foregrounds can be perfectly subtracted on 
scales k > kmm- However, even on these scales, temperature fluctuations due to Galactic syn- 
chroton, extra-galactic point sources and free-free emissions can still significantly overwhelm 
the 21cm signal, and its subtraction can remove powers that carry essential information from 
primordial non-Gaussianity. Foreground cleaning by polynomial fitting (see e.g. [43, 44] 
and references therein) could also possibly "fit out" the non-Gaussian features, and a non- 
parametric cleaning scheme [45, 46] may be a more appealing alternative. It is hoped that 
by the time the SKA comes into operation, its precursors [47-49] would have shed enough 
light on the nature of foregrounds that they may be feasibly subtracted. 

6 Conclusions 

We have shown how information on primordial non-Gaussianity can be extracted from the 
21cm radiation signals - a connection which deserves to be explored more deeply. Our primary 
focus was to study the manifestation of various shapes of non-Gaussianity in the 21cm power 
spectrum, and thus evaluating the prospects for future radio telescopes in detecting these 
signatures. Our calculations are based on recent analytic work on the bias of large-scale 
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structures, modified and extended to scales where there are good prospects of detecting the 
21cm signals from the epoch of reionization. We found that the shapes of non-Gaussianity 
manifest as different scale dependence in the bias of ionized fraction as shown in Fig. 1. 

Our analysis demonstrated that the /nl constraints from the SKA will be competitive 
with those from upcoming CMB experiments such as Planck (and much better than current 
LSS constraints), whilst futuristic radio arrays can tighten the /nl constraints down by 
at least an order of magnitude towards o"(/nl) = C'(l). There are several ways in which 
the constraints can be further improved. Firstly, we note that our analysis uses only the 
information in the power spectrum of the 21cm radiation, whereas much stronger constraints 
are expected from measurements of the bispectrum [9]. Secondly, the 21cm measurements can 
be made across a larger number of redshift bins. Indeed, 21cm tomography has been shown 
to potentially outperform CMB experiments (which provide measurements of fluctuations at 
a single redshift) in constraining key cosmological parameters [35]. Finally, the CMB, large- 
scale structures and 21cm signals probe non-Gaussianity on different physical and time scales, 
and combining these constraints will be crucial in the understanding of complex physics such 
as scale-dependence and time-evolution of /nl- 

A detection of non-Gaussianity will clearly have far-reaching consequence for our under- 
standing of cosmology, and the power of the 21cm signal in probing non-Gaussianity is only 
becoming better understood within recent years. With future work, we plan to show how the 
21cm radiation can be used not only to constrain but also to disentangle various shapes of 
non-Gaussianity. This will be crucial in the understanding of the physics underlying inflation. 
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A 21cm power spectrum with the Sheth-Tormen formahsm 

We outline a method of extending our calculation of the ionized- fraction bias b^, based on 
the Press-Schechter collapse formalism, to the Sheth-Tormen formalism. The idea is based 
on averaging the large-scale bias over the mass function: 

bsT = 7oo , , , , (A.l) 

Note that these equations reduce to the PS case if p = and q = 1- 

To calculate bx in our chosen reionization model (3.7), we also need the collapse fraction, 
/coll, which, in the Sheth-Tormen formalism is quite complex (see the expressions derived 
from excursion-set theory in e.g. [50, 51]). However, the scaling 6c — )• ^/qSc was found to be 
a reasonable approximation for ellipsoidal collapse integral to the Sheth-Tormen approach 
[52, 53]. 

From [29], the bias shift is given by 
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A6sT(fe,M) 



+ 



2 1/2(1 + (gi/2)p) 



I{k, M) 



q p 

2 ^ 1/2(1 + (gi/2)P) 



(A.4) 



din (Tm 



These expressions together with the Alvarez bias (4.18) yield the ionized-fraction bias, b^. 
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